(N 



cd 



(N 

0^ 



(N 



■ r^)'.-' ' ■ -ti Availableonlineatwww.sciencedirect.com . 

y^"^ ••;* ScienceDirect AnnaisofPhysks 



ELSEVIER Annals of Physics 00 (2013) l- fTe] 



A quantum dot close to Stoner instability: the role of Berry's Phase 

ARIJIT SAHA AND YUVAL GEFEN 

' _ ' Department of Condensed Matter Physics, Weizmann Institute of Science, Rehovot 76100, Israel 

CN IGOR BURMISTROV 

> 

^ Landau Institute for Theoretical Physics, 1 19334 Moscow, Russia 

^ ALEXANDER SHNIRMAN 

C^ Institutfiir Theorie der Kondensierten Materie, Karlsruhe Institute of Technology, 76128 Karlsruhe, Germany 

^^^ 6- DFG Center for Functional Nanostructures (CFN), Karlsruhe Institute of Technology, 76128 Karlsruhe, Germany 

^ ALEXANDER ALTLAND 

i-i^ Institutfiir Theoretische Physik, Universitdt zu Koln, D-50973 Koln, Germany 



Abstract 

The physics of a quantum dot with electron-electron interactions is well captured by the so called "Universal 



I 

r^ Hamiltonian" if the dimensionless conductance of the dot is much higher than unity. Within this scheme interac- 

O tions are represented by three spatially independent terms which describe the charging energy, the spin-exchange and 

^ O^ the interaction in the Cooper channel. In this paper we concentrate on the exchange interaction and generalize the 

functional bosonization formalism developed earlier for the charging energy. This turned out to be challenging as 
PsJ the effective bosonic action is formulated in terms of a vector field and is non-abelian due to the non-commutativity 

¥^ of the spin operators. Here we develop a geometric approach which is particularly useful in the mesoscopic Stoner 

regime, i.e., when the strong exchange interaction renders the system close the the Stoner instability. We show that it 
is sufficient to sum over the adiabatic paths of the bosonic vector field and, for these paths, the crucial role is played 
■^1- by the Berry phase. Using these results we were able to calculate the magnetic susceptibility of the dot. The latter, 

^^ in close viciiuty of the Stoner instability point, matches very well with the exact solution (Pis'ma v ZhETF 92, 202 

O (2010)). 



Keywords: Quantum Dot, Berry Phase 



S^ 1. Introduction 

Cu Over the past few decades physics of quantum dots (QDs) has become a focal point of research in nanoelectronics. 

The introduction of the "Universal Hamiltonian" fll El |3] ID has made it possible to take into account the effects of 
electron-electron (e-e) interaction within a quantum dot (QD) in a controlled way. This approach is applicable for 
a normal-metal QD in the metallic regime when the Thouless energy Eji, and the mean single particle level spacing 
6 satisfy g = Erh/S » 1 (^ is the dimensionless conductance). Within this scheme interactions are split into a sum 
of three spatially independent contributions in the charging, spin-exchange, and Cooper channels. The charging term 
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Figure 1: (Color online) Cartoon scheme of an isolated QD . Here d is the single particle level spacing, (Do is the zero component of the auxiliai'y 
HS vector bosonic field, and B is the applied magnetic field. 



leads to the phenomenon of Coulomb blockade, while the spin-exchange term can drive the system towards the Stoner 
instability |5|. 

In bulk systems the exchange interaction competes with the kinetic energy leading to Stoner instability. In finite 
size systems mesoscopic Stoner regime may be a precursor of bulk thermodynamic Stoner Instability IJJ. More 
precisely, one distinguishes three regimes depending on the strength of the exchange interaction: (a) a phase with 
the total spin of the dot equal zero, (b) the mesoscopic Stoner regime in which the total spin of the dot is finite but 
not proportional to the volume of the dot, and (c) the thermodynamic ferromagnetic phase where magnetization is 
proportional to the volume. The mesoscopic Stoner regime can be realized in QDs made of materials close to the 
thermodynamic Stoner instability, e.g., Co impurities in Pd or Pt host, Fe dissolved in various transition metal alloys, 
Ni impurities in Pd host, and Co in Fe grains, as well as new nearly ferromagnetic rare earth materials l'6]|7][8l. 

Notably, the inclusion of the spin-exchange turned out to be non-trivial as the resulting path integral action is 
non-Abelian |9 10 1. To understand the complexity of the problem we compare with the case when only the charging 
interaction is taken into account 1,1 li fT2l [T3l . It was suggested by Kamenev and Gefen | ,11J to take the following 
steps in solving that problem: (a) start from a fermionic action which includes an e-e interaction term quartic in the 
fermionic Grassman variables, (b) perform a Hubbard-Stratonovich (HS) transformation by introducing an auxiliary 
bosonic field, (c) perform a gauge transformation over the Grassman variables which makes all the non-zero Matsubara 
components of the (HS) field decouple from the fermionic fields in the action, and, finally, (d) integrate out the 
fermions. The resulting, purely bosonic action, is quadratic in the bosonic non-zero Matsubara components, which 
renders the problem easily solvable. The trick of gauge-integrating over Grassman variables does not work for the 
non-abelian case |]9l, so that an alternative approach is needed. 

There have been several attempts to account for charge and spin interactions in a QD including a rate equation 
analysis 1.14.. 15J and a perturb ative expansion [9J. AUiassid and Rupp 1 14| have analyzed some aspects of the problem 
exactly. More recently an exact solution of the isotropic spin interaction model based on the generalized Wei-Norman- 
Kolokolov method |16| has been presented |10|. In this exact solution several observables, including the tunneling 
density of states and magnetic susceptibility have been calculated below the Stoner instability point for an equidistant 
spectrum. The effects of disorder have been addressed in Ref IfTTl . The tunneling density of states exhibits a non- 
monotonous behavior as a function of energy, and the magnetic susceptibility emerges out to be a sum of Pauli and 
Curie like terms. 

In this article we present an approximative geometric approach to tackle the isotropic spin-exchange model. Our 
results are in agreement with the exact results IfTOl for the partition function and the magnetic susceptibility within 
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the mesoscopic Stoner instability regime. Our rationale behind developing an approximation scheme, given the exact 
solution, is the high complexity and inflexibility of the exact method. We thus expect to be able to apply our geometric 
approach in cases where the exact method is inapplicable or too complicated. These should cover a broad range of 
problems involving spin transport through QD coupled to normal leads and through an array of QDs . 

This paper is organized as follows. In Seep] we consider an isolated QD with isotropic exchange interaction, 
and derive an effective action in terms of an auxiliary HS bosonic vector field. In SeclSlwe perform a perturbative 
expansion of our effective action in powers of the Berry's connection operator, and show that in the lowest order 
of this expansion the Berry phase governs the effective action. In Sec|4]we calculate the partition function and the 
magnetic susceptibility of the QD , and discuss the effect of Berry phase on susceptibihty. SeclSlcontains a summary 
of our analysis. In the appendices (A, B and C) we present alternative methods of calculation, which provide further 
justification of our results. 

2. Hamiltonian and Effective action 

In this section we perform a HS transformation and obtain the effective action in terms of an auxiliary HS bosonic 
vector field <P. We, then perform a unitary rotation R to the instantaneous direction of $ and rewrite the effective 
action in terms of the Berry connection operator R^^R. Finally, integrating out the fermions we obtain the effective 
action of the isolated QD in terms of the zero component of the HS field and R^^R. 

A quantum dot in the metallic regime, g » 1, is described by the universal Hamiltonian |[T1: 

H^Ho + Hc + Hj+ H,. (1) 

The noninteracting part of the universal Hamiltonian reads 

Ho = J] eaal^a^^^ , (2) 

where e^ denotes the energy of a spin-degenerate (index cr) single particle level a. The charging interaction term 

Hc = Ec{N-Nof (3) 

accounts for the Coulomb blockade. Here, Ec denotes the charging energy of the QD, A^o represents the background 
charge, and N - Yja.a '^I.o-flff.o- is the operator of the total number of electrons of the dot. For the isolated QD the total 
number of electrons is fixed and, therefore, the charging interaction term can be omitted. The term 

Hj = -/S' (4) 

represents the ferromagnetic (J > 0) exchange interaction within the dot where S - Yja "a,a-,Sa-,o-,aa,,o-2 is the operator 
of the total spin of the dot. Here 8(^,0-2 = (l/2)(?o-icn, where cr = (o-x,o-y,o-^) is a vector made of Pauli matrices. The 
interaction in the Cooper channel is described by 

a 

In what follows we do not take into account i/^ for the following reasons. For the dots fabricated in 2D electron gas 
the interaction in the Cooper channel is typically repulsive and, therefore, renormalizes to zero |2 1. In the case of 3D 
quantum dots realized as small metallic grains, the interaction in the Cooper channel can be attractive, giving rise to 
interesting competition between superconductivity and ferromagnetism lfT8l[T9ll20ll . In that case we assume that there 
is a weak magnetic field which suppresses the Cooper channel. 

As explained above we restrict ourselves to a simplified version of the universal Hamiltonian, where the interaction 
in the charging and Cooper channel is set to zero (more precisely, the charging energy is fixed because the number of 
particles is fixed): 

// = ^ e^al^^a^,^ - JS^ . (6) 



by 
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Our aim is to calculate the partition function Z, - J D^D^ exp [»S>p], where the imaginary time action is given 



p p 



P 



dr 



2_^ ^aa-i-dr - £a + ^l)^aa- + J /_J '/'«<r, So-.o-,'/'^ 



(Ti (T2 V' O-CTt 



(7) 



Here fj. is the chemical potential, /3 = l/T, T the temperature, and we have introduced the Grassmann variables 
^ff = (lAffT' ^a[f , ^a - (^oi, 4'al) to represent electrons on the QD. 

In Eq|7] the exchange energy is quartic in the fermionic fields. Hence, we can perform a HS transforma- 
tion to obtain an effective action quadratic in the fermionic fields, with an auxiliary vector bosonic field 0(r) = 
(<1)j;(t), <1>v(t), ^^(t)). The effective action reads 



5^,o= ffl'T Y^^>„[-^,-6,,+^i-6■§)^',-^^ 



(8) 



Integrating out the fermions, we obtain the effective action in terms of the auxiliary vector field <P only 



vn / ^ ^. r |(D|2 

S^^Y^Xx \n[-dr -6„ +//-(!)■ S) -J dT^-^ 
" 

p 

= ^tr lnl-5^ -e„+fL- (I)(t)«(t) ■ -J - \ dr 



02(t) 
47 



(9) 



Here n{T) is a unit vector along the direction of <P(t) and <1)(t) = \^{t)\. The first part of the action in Eql9^ describes 
the coupling of non-interacting electrons to a time varying magnetic field (exchange field) of magnitude <t)(T) and 
direction «(r). The resulting bosonic action in Eq|9]is non-abelian due to the non commutativity of the Pauli matrices. 
Note that the problem is 'isotropic' in the sense that the outcome should not depend on the initial and final direction 
of ft. 

We are guided by the idea that close to Stoner instability the amplitude of the exchange field <I> is large, i.e., a 
large total spin develops on the dot (<P ~ 7S). In that situation one can distinguish between the adiabatic and the 
non-adiabatic paths of *5(t). The former involve Matsubara frequencies such that \(jl>„,\ <k (D. We argue that the non- 
adiabatic paths do not contribute considerably to the partition function (except for providing for proper normalization), 
since the electrons do not manage to react to the fast changes of iP. Hence, we concentrate on the adiabatic paths and 
perform an expansion in the time variation of n. We transform to a coordinate system in which it coincides with the 
z-axis 



Ra-,Rr 



(10) 



where 7? is a unitary rotation matrix. Eq. ( 10 1 identifies R as an element of SU(2)/U(1). Indeed, if we employ the 
Euler angle representation 



R - exp 



r id, 1 




r io 1 




r ill, ] 


[-T^\ 


exp 


[-2"-^ 


exp 


[-y'^^J 



(11) 



then the angles and 6 determine the direction of «, while ij, is arbitrary, i.e., the condition ( 10 1 is achieved with any 
value of if,. Thus, ij, represents the gauge freedom of the problem. We obtain 

Y^ tr In i-dr -e„+ii- <^{T)ft{T) ' 1 1 = J] ^'" '" ("^^ -£a+^i- *&(t) y - R'^drR^ . (12) 
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In the transformation ( 12 1 one can, in the spirit of Ref. ifTTI . think of R as being appHed to the fermionic field. That is, 
one first introduces ^''^ via ^'a = R^'[^ and, then, integrates over ^''^. In this case it is convenient to choose the gauge 
ifr so that R remains periodic in Matsubara time upon a continuous change of (pto (p + In. Then ^^ is anti-periodic as 
it should be. This can be achieved, e.g., by fixing the gauge as i//{T) - -(pir). In what follows we work in this gauge. 
Next, we represent the amplitude of the exchange field as a sum of its zero frequency component and of the 

rest, ^(t) = 00 + d^ir), such that J S<t>(T)dT - 0. It is clear and can be easily checked that adiabatic longitudinal 



fluctuations (50 do not contribute substantially to the effective action (see Appendix A for a more formal treatment). 
Thus, we disregard that part of (50 and obtain 



2 tr In (-5, - e„ + /i - (1)(t) ^ - R'drRJ * ^ tr In (-5, - e„ + /z - Oq ^ - 77(t)) 



(13) 



where 



1](t) = R~'(drR) 



I ■ I ■ T 1 ' ■ r 

— —-(f>(cos0- l)cr. + - sin0(cos0crv- + sin^cr^) - t 6 



cos mcTy - sm 



bo-.,] 



is the Berry's connection operator in the gauge ipir) = -0(t). 



(14) 



3. Expansion in powers of the Berry connection R ^R 

In this section we perform a perturbation expansion in powers of the operator rj, and obtain an effective action which 
encompasses both the longitudinal fluctuations and the Berry phase term. 



We aim at the expansion of the action d9]l in powers of the Berry connection operator rj - R ^R (see 14 1. It is 
expected that this expansion quickly converges for the adiabatic paths of 0(t). We write iSo = Zi„>S^'\ The zeroth 
order term S)^^ can be obtained by calculating the grand canonical potential of the noninteracting electrons subject to 
a constant Zeeman field Oq- We obtain 



S%^ ^Yj^r\n i-dr - e„ + ;/ - Oo ^) - { dr ^ = -y6Qo(<l>o) " T7 " Z 5 ^'^'«^*^-'« ' ^^^^ 



where 



-/JQo(<l>o) = InZo = 2 [ln(l + g-^'i^"-^-/')) + in(l + ^-Z^l^"^^-")) 



(16) 



Here Zq is the partition function of noninteracting electrons subject to a magnetic field of amplitude Oq. To determine 
the Oo-dependent part of Qo we calculate 



r(<Do) = 2 



50 



^?I4 



1>o\ .L Oo 



(17) 



where ^„ - ^a - 1^ and /(e) = (expljSe] + 1)"' is the Fermi distribution function. At zero temperature r(<l)()) is the 
number of orbital levels between n - 00/2 and ^i + <l>o/2. Assuming a constant density of states (equidistant spectrum) 
we obtain 

5Qo 1 - - v02 



r = vO, 



'0 : 



vOn , Qo — const. ■ 
50o 2 4 



(18) 
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Strictly speaking (18 1 is valid at temperatures higher than the level spacing 6, i.e., for T » (5 = v . At lower 



temperatures step-like dependencies are expected. Yet, in a "coarse-grained" sense (18 1 holds at lower temperatures 
as well. Finally, 



P^l sr P 



S'^^^ const. -"-^-Y.f-j^^nS^^-.n, 



47* ^ 47 



(19) 



where 1 /7* = 1 /7 - v. As we consider the regime close to Stoner instability, we have J* :s> J > 0. Realistically, 
the quantum dots are disordered and the assumption of an equidistant spectrum is too naive 1.171 . Due to disorder we 
should have 

Qo = const. - —^ + 6Slo , (20) 

where v is the average density of states. This question was originally addressed by Kurland et al. 1 1 j and was recently 
analyzed by Burmistrov et al. IfTTl . Roughly, 6£lo is of order ±(9(l)Oo. In the present paper we disregard disorder. 



3.1. The first order contribution S^ 



In the first order in rj (Eq. 14 1 we obtain 



^o' = 'l j dTtr[G\T,TUT)] = -^tr G"(e„)^ J drTjir) 



(21) 



where G^^,(e„) = 6a.a' [isn -^ar - y Oq) andfi„ - 7r(2n-i-l)/j8 are the fermionicMatsubara frequencies. Calculating 
the sum over b„ we obtain 



^:=4Z" 



/i..4)./k-^)-f„4/L.^)-/k-^ 



\\cr, I dT rjij) 




.(22) 



From tr(;7(T)) = we conclude that 



^ Jo -^ Jo 



(23) 



where F was defined in Eq. (17 1. The contribution to the effective action, given by Eq 23 is proportional to the Berry 



phase. The coeflicient in front of the Berry phase F(<t)o) is roughly equal to the the number of single occupied levels, 
i.e., the number of uncompensated spins which acquire the Berry phase. 

3.2. The second order contribution S}^ 

To calculate the second order contribution we introduce the notations G?,, = (ie„ - ^a + 'l*o/2)"' and cr-^n = 
1(1 + cr-,) and we obtain 

^''^^^-\YjYj ^' ('?'«^«+'«'?-'"^«) ^-\YjYj^' (^m(G?crT + Gjcri)„+^/7-m(G?crT + Gjcr^),,) . (24) 

a n.m a n,m 

Here rj^ = \ J dr ri(T)e"^"'^ and oj,,, - Imn/jS are bosonic Matsubara frequencies. For the transverse, i.e., the spin 



flipping part of //, we obtain 



,(2) ^ ^ V 



■^O.i 



2 ^ 



tr('7mCrt7;_mCr|) 



1 



On - iaj„ 



+ tl"('?mCr|77_,„Cr|) 



1 



On + ico,, 



(25) 
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In the adiabatic limit |w,„| <K Oq this gives 

■^f± ~ :ir y tr(?7mcr|7;_,„cr|) = — I dr ti {r](T)o-^T](T)a-i) . (26) 

•Po ^ ^0 J 

Hence, substitution of Eq{T4]leads to the result 

\ drW + (sme4>f\^ {dTif. (27) 

J ^ J 4<l)o J 



^(2) ^ _ 
t.-L 4(Do 





where n is the angular velocity describing the motion of the magnetization vector O on the sphere. 
For the longitudinal part of rf we obtain 

(2) ySv V I r^ . ^ 

^V\\ ^ ~2 [^i'('?'"=o'^T^"'=oO"t) + tr(?7m=oo-i'7m=oo-|)] = -— II dT(f>{\- cos0)l . (28) 



Finally ^S^' = »S^ j^ + S)^... Both this terms are ~ liv{nf-. In comparison S)^^ of (23 i is ~ pv(t>i)n. Thus for all adiabatic 
frequencies |w,„| <s Oq the Berry's phase action S^^ dominates. In what follows we restrict ourselves to the first order, 
Berry phase, contribution S]^ . 

4. Partition function and magnetic susceptibility 

In this section, after performing the path integration over the adiabatic paths ofn(T), we obtain the partition function 
of the problem as a function o/Oq only. This allows us to calculate the magnetic susceptibility consisting of the Pauli 
and Curie terms. 

Above we have obtained the following effective action for the adiabatic paths 

S^ ^ const, -j^'l^f-j S^m5^-m +- \ dT(\-cos0)(l>. (29) 

This action governs fluctuations at low frequencies |(w,„| < Oq. Thus, at these frequencies the system reduces to a large 
spin of amplitude ~ F. Indeed the last term of the action (J29]l is just the well known Wess-Zumino action of a free spin. 
While a partition function of a true free spin is trivial to calculate, our spin "lives" only at the adiabatic frequencies. 
We, thus, define our functional integral for the partition function as an integral over all paths whose frequency scale 
is cut off by Oq. We have 

' " ' (30) 



Z = N f D^<^ e^" , where £)^(D = ]~[ c/^O, 



In Eq 30 the index / defines a lattice partition of the interval [0,/?] defined so as to limit frequencies to values \lo,„\ < Oq. 
The measure c/^O,- is an ordinary Cartesian measure. The Hubbard-Stratonovich transformation and integrating out 
the fermions produce the following normalization factor 

N^n{^) (31) 

Here N - l/\A'^^2)4?{'^'^2)\ and A^ - %/27iT » 1 is the number of positive Matsubara frequencies taken into 
account (those that are below the cut off Oo). 
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4.1. Calculating the path integral 

We are guided by the idea that the most important paths are those of almost constant radius <1>(t) « Oq. Indeed, 
as can be seen from ( 29 1 the longitudinal fluctuations dO are suppressed by a factor containing the bare J. By 
contrast, close to the Stoner instability, the zero frequency amplitude Oq is only weakly suppressed {J* » J). The 
larger is Oq, the bigger is the phase volume of possible transverse fluctuations. The latter are only "penalized" by 
the Berry phase term. An entropic "phase space" argument, outlined above, reveals that (Dq should assume a large 
value. Below we find out that the typical value of Oq is of the order of 7*, whereas we have |(JO„,| ~ VtV and 



|5<I)(r)| ~ V 2m l'5'l'ml ~ yJ*J ^ Oq. To derive all these, we begin by converting the measure to a polar one, which 
is then adjusted to an integration over paths of almost constant radius Oq: 



n 



©3,5 ^ I I (t,fd(bidiii 

ietime intervals 



^2Z.i„a.,j-[^^,^._^. 



2(27V+1) -i2 2,'50T 
^ ° 



ri d^idfii 



(32) 



where in the last identity we have dropped the term in the exponent 






O^A 



p 
iA^J 



dT(6^(T)f ^ 



\(jJ„i\<lL)i 



Yj \^^" 



(33) 



Here w^. ~ 1 /At is the ultra-violet cutoff". Taking consistently the adiabatic cutoff o), 
is negligible in comparison with the 



i&o we observe that this term 



PI J term in the action (29 1. 

We now perform the Gaussian integration over the longitudinal fluctuations 5<I>„„ the net eff'ect of which is the 
partial cancellation of the normalization factor (31 1 by (j3/J)^. We are left with 



z 



^*(ir(f)7*>' 



J 



exp 



47* 



(D;: 



p 

r ^ iT r . 

I Dn exp — I iiT0(l - cos^) 



(34) 



We next proceed to the integral over the transverse fluctuations, i.e., the fluctuations of n(T). These fluctuations are 
"penalized" by the Berry phase term in the action. The Berry phase term is given by the solid angle swept by the n{T) 
path. Since F ~ vOq s> 1 (this has to be checked for self consistently), these solid angles must remain small, and we 
can restrict ourself to a Gaussian expansion around a static value of n, namely no, the relevant paths will be given by 
the small variations around this static «o. At the end we should average over all possible directions of «(). 

The Gaussian integration is easily performed if we note that the integration measure D« is given by H/ d4>id{\ — 
cos Oi). Thus, introducing y =\ - cos 9 we obtain 

^m exp I [ dT4>{\ - COS0) * ^DcjiDy exp ^ f ofr^y * Y[ (-^\ ■ (35) 

L J L J "'=^ 

In evaluation of the integral we extended the integration limits of both and y to [-00,00] even though, e.g.,3; e [0,2]. 
This is justified for F » 1 . Thus we obtain 



z 



*(irr7)7*Mfr'-*^n(,-i:)'=*(fr(?)U)7*>-*-^n 






a)„ 



(36) 



The product appearing in Eq. 36 is cutoff at m - N, i.e., at ojm - ^o- This hard cutoff" is an artifact of our rather 



hand-waving approach. Within this method we have no way to determine what kind of cut off' should be employed. 
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One possibility, which is supported by the calculation in Cartesian coordinates provided in [Appendix B 
soft cut off. This gives 



IS to use a 



" '»»>'- 



nehnf^" 



m=l 



m=I 



2 \ 



sinhf^-f^l expf^-f^l 



m) 



(37) 



Evidently, employing a different cut off, one would obtain a result which is different in the exponential. For example, 
the hard cut off gives ^^ in the exponent. This would slightly modify the numerical coefficients in the final result 



for the susceptibility. We keep here the soft cut off result, as it is supported by the calculation in Appendix B and 
provides an excellent approximation to the exact solution 1 10|. 



We approximate the ( 1 /vJ)'^ factor in ( 36 1 by 1 , which is legitimate near the Stoner transition. Finally, integrating 
over (Do we obtain 



Z^N 



(^"i5ri 



ii<l)o4;riI)o exp 



47* 



sinh[^] 



2 



4jiN 



(i) 



N/j* 

7 



3/2 



exp 



jS7* 



(38) 



4.2. Magnetic susceptibility 



We obtain the magnetic susceptibility using Eq 38 as follows 



X 3 



1 dlnZ 1 



y6 



1 



dJ 2(1- Jv) 12 (1 - v7)2 " 
We observe that it consists of a Pauli and a Curie contribution. In comparison the exact solution of Ref. ITOl reads 



(39) 



J* 



Zexact - I y I exp 



3/2 



fi{r - J) 



and the magnetic susceptibility 



Xexact 



2(1 



- Jv) 12 



1 



(i-vy)2 



(40) 



(41) 



J the extra -1 factor mxexaa is immaterial and we obtain an 



In close vicinity to the Stoner instabihty, i.e., 6 ~ v 
extremely good approximation to the exact result. 

The Pauli-like (with an upward renormalized g-factor) susceptibility (first term in Eq. (39l) dominates when 
T » J*. In the low temperature regime, T <K J*, the Curie-like part (second term) dominates. In this regime 
the average spins scales as V< S >^ ~ J* 15. This Curie like contribution in the magnetic response can be tested in 
materials close to the Stoner instability such as Pd (JjS - 0.83) or YFe2Zn2o {J 1 6 - 0.94) |8|. It is important to 
understand that the Curie part of the susceptibility represents a mesoscopic effect. The density of states of a QD v 
scales linearly with the volume of the dot, v ~ V. On the other hand J ~ l/V. Hence the Pauli like part of the 
magnetic susceptibility is proportional to the volume V and is an extensive quantity. On the other hand, the Curie 
susceptibility is intensive as vJ is scale invariant. Therefore for a fixed temperature T, if one gradually increases 
the size of the system, the Pauli part grows and, eventually, the Curie susceptibility becomes negligible compared to 
the Pauli one. This behavior is shown in Fig|2ta) where vJ is chosen to be in close vicinity of the Stoner instability 
regime. To illustrate the nature of Pauli and Curie like suceptibility more, we show their behavior as a function of 
the Stoner parameter vJ in Fig.|2|b). It is clear that for a fixed dimensionless parameter y37, the Curie susceptibility 
(green curve) becomes more dominant over the Pauli one (red curve) as we approach towards the Stoner instability 
point. On the other hand, the blue curve shows the behavior of Jxexact which in close vicinity to the Stoner instability 
point matches quite well to our approximate answear (Eq. ( 39 1), shown by the magenta curve in Fig. l2lb). 
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Figure 2: (Color online) (a) Magnetic susceptibility of an isolated QD in the T - V plane where T is the temperature and V is the volume of the 
background 3D material. Here v] = 0.9 i.e. close to the Stoner instability point, (b) Dimensionless susceptibility Jx is shown as a function of the 
Stoner parameter vJ. Here the dimensionless parameter /?/ = 1 .0. The red, green, magenta and blue curve correspond to the Pauli like susceptibility 
(first term in Eq. |39}), Curie like susceptibility (second term in Eq. J39l), our approximate expression for the magnetic susceptibility (Eq. J39}) 
and the exact expression of susceptibility (Eq. J4T)) respectively. 



5. Summary and Discussion 

In this paper we have considered an isolated QD with an isotropic exchange interaction. In the path integral 
HS formulation this problem leads to an effective non-abelian action. Here we have presented an approximative 
geometric approach in which the Berry phase controls the dynamics of the direction « of the HS magnetization vector 
O. Close to the Stoner limit, i.e., for v7 — > 1, our approach reproduces well the exact solution of Ref. IfTOl which 
comes at the expense of a hard calculation based on the generalized Wei-Norman-Kolokolov method fTSI. Note 
that even zero frequency observables (for e.g. susceptibility), involve summation over finite frequency fluctuations 



(Eq. (B.14i), including both low and high frequency contributions. Although the exact result describes both the low 
and high energy contents of susceptibility, we here show that if one is interested in the low energy regime of long 
range fluctuations, the information of the exact result can be obtained from our physically motivated and user friendly 
approach of the invariant action. However if one insists on knowing ultraviolett behavior, then one may resort to an 
equally simple gaussian expansion around stationary points. Hence, our approximative geometric approach covers 
much of the contents of the observable, in a manner suitable for further generalization. Therefore, we believe that 
this approach could be very useful in situations in which the exact method is inapplicable. These may be, e.g., the 
problems of charge and spin transport via quantum dots coupled to normal or ferromagnetic leads. 

Strictly speaking our results are valid for <5 = v' ~ J < T < J* . Yet, for the equidistant spectrum assumed in this 
paper we do not expect major changes at lower temperatures. In the presence of disorder the low temperature regime 
T < J* is itself more subtle (TJ]. 



Finally, it is important to mention that, as discussed around Eq. ( 20 1, disorder can influence the grand canonical 
potential of non-interacting electrons Oq quite essentially. Yet, it is easy to show that the Berry phase part of the action 
is quite insensitive to disorder 



6. Acknowledgments 

This work was supported by GIF, Einstein Minerva Center, Sonderforschungsbereich TR 12 of the Deutsche 
Forschungsgemeinschaft, EU FP7 grant GEOMDISS, the Russian-Israel scientific research cooperation (RFBR Grant 



Saha, Gefen, Burmistrov, Shnirman and Altland / Annals of Physics 00 (2013) 1-\16\ 



11 



No. 1 1 -02-92470 and IMOST 3-8364), the Council for Grant of the President of Russian Federation (Grant No. MK- 
296.2011.2), RAS Programs "Quantum Physics of Condensed Matter" and "Fundamentals of nanotechnology and 
nanomaterials", the Russian Ministry of Education and Science under contract No. P926. 

We acknowledge useful discussions with Gabriele Campagnano, Igor Lerner, Mikhail Kiselev, Jiirgen Konig and 
Alexander Mirlin. We are grateful to Ganpathy Murthy for providing us with notes of his calculations on "Universal 
interacting crossover regime in two-dimensional quantum dots" (Ref. l.21J ) and a detailed explanation. 



Appendix A. Unimportance of longitudinal fluctuations 

Here we show why the longitudinal fluctuations of 6^ can be disregarded. 

Following the spirit of Ref. ifTTI we gauge out in Eq. (12 1 the fluctuations 50 by using a (non-unitary) transforma- 



tion V - exp -f ''^z ' where ^(t) = J 6€>(T')dT'. We obtain 



J^ tr In {-dr -ea+H- <1>(t) ^ - R'drR^ = X ^'" ^" [~^' - 6„ + // - Oq ^ - //(t)) , 



(A.l) 



where 



rj{r) = V-'R-'(drR)V 



I . I 

•-(^(cos 6 - l)cr. H — exp 



(0[sin0( 



cos (per, + sm (Act, 



,)]-«[ 



cos (pcTy - sm 00" V ) exp 



-1) 



o:, 

2 " 



(A.2) 



This means that, in fact, //(r) appearing in (14 1 is the expression (A.2 1. We observe that the longitudinal part of 77 of 
(A.2 1, which is responsible for the Berry phase term in the action, does not contain the operators V, V ' and, thus, 
is not influenced by the longitudinal fluctuations 50 which appear only in the factor ^. Moreover, substituting 77 of 



( A.2 1 into Eq. ( 26 1 we observe that in the adiabatic limit the longitudinal fluctuations drop out also in the second order 
terms. Thus, disregarding the longitudinal fluctuations in (fT3]l, and the factor V altogether, was justified. 



Appendix B. Expansion in smaU transverse fluctuations 

Here we present an alternative method of calculating the path integral in Cartesian coordinates. The advantage is 
the higher level of accuracy in handling the integration measure. An adiabatic cut off need not be postulated here. 
Rather, a soft adiabatic cut off appears naturally. 

Appendix B.l. Effective action 

We start from the Hamiltonian (|6]l and rewrite the action (|7]) in the Matsubara representation as 



S^v^fi 



2_] i/aa^niif^n ' ^a + lj)^ao-,n + J /_j 



7 . Yao-x^o-yiTyYao-y 



7 . Yao'^^^o'xO'^.Yaa'y 



(B.l) 



We apply the HS transformation on Eq B.l to obtain an effective action quadratic in the fermionic fields, with an 
auxiUary vector bosonic field (O) for the spin degrees of freedom. Hence, the effective action reads 



Oiy^O — p 



2:^aa;n(f^n ~ ^a + lj)^aa;n ~ /i "^ ' /, '^a(T,^a-,tT2'l^aa-2 ~/. 



^,A-„ 



4/ 



where the Matsubara expansion for the bosonic HS real vector field O reads 



(B.2) 



(B.3) 



m^O 



Saha, Gefen, Burmistrov, Shnirman and Altland / Annals of Physics 00 (2013) 1-\16\ 



12 




\50, 



X 



Figure B.3: (Color online) Schematic of the auxiliary bosonic HS field <5 where <Po is the zero-component and (SO,, (5<I)y,i5<D; are the non-zero 
Matsubara components. Among the three non-zero Matsubara components S(^- is the longitudinal and (JO,, (SO,, are the transverse fluctuations 
respectively. Here in our analysis the longitudinal component 6^. is chosen in accordance with the direction of <5o and S^x, S^v are the two 
transverse componets. 



Here Oq is the zero component of the bosonic HS real vector field <5 and (JO^ = O^ are the non-zero ones. 

Our strategy is to, first, integrate over the non-zero Matsubara components (50,„ while keeping Oq fixed and, then, 
integrate over all possible <Po. The vectors (JO,„ have three components among which one is longitudin al an d the other 
two are the transverse with respect to the current direction of <Po, which is schematically shown in Fig 
we choose the basis for (5<5„, in accordance with the direction of Oq, i.e., the axis z for 5<5,„ is along Oq. 



B.3 



Moreover 
In other words 

for the time being we break the symmetry of the isotropic problem choosing a particular direction of Oq and select 
the basis for (50,„ in accordance with the current direction of Oq- Note that at the end of the day one should integrate 
over all <Po in order to restore the global symmetry of the problem and obtain an expression for the partition function 
and the susceptibility compatible with Eqs|38]and 39 In terms of <^o and the cartesian fluctuations around it, Eq B.2 
can be rewritten as 



kSiy^o — /? 



02 



Tropin 2_j ^a-m [(Go,L,„,)^«i,«2 - ^•^(«l-«2) ■ S] ^l/a,n2 ' Tj ~ Z_i 



5<b,„6<i>-„ 



mi:!.) 



47 



(B.4) 



where G^^^^ - is„ - {ea, - lu) - (^oS z is the single particle Green's function of the electrons subject to a constant 
magnetic field (Sq. For the partition function we obtain 



1 \^('^+3)/«\3/2/flvW r 



Integrating out the fermions we obtain 



Z(m)^N 



-.rlPVUH 



6)(?)7 



D^e^" , 



where the effective action >So depends on Oq and 6^m and can be written as 






(B.5) 



(B.6) 



(B.7) 



m#0 
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Here tr stands for combined time (Matsubara frequencies) and spin trace. 

Appendix B.2. Expansion to the second order in fluctuations 

We expand the effective action (Eq B.7i up to the second order in the fluctuations 5<5„,. The zeroth order con- 
tribution has already been calculated in Sec. [3] It is easy to show that the first order contribution vanishes and we 
obtain 



S, . const. - ^ - y ^-^^^ - X ^ tr fG„(^<5 ■ S)G,(66 ■ S)] 



47' 



into 



After a straightforward calculation this gives 



^ ^ tr [Go(6^ ■ S)Go(66 ■ S)] = /3r ^ 



47 



50_.„5<D+._„ 



(B.8) 



m#0 



Oo - /<y„ 



2 



O, 



pr 






'i>{) 



>0 



(B.9) 



where 5<^± - (iSO^ + /(J<l)y)/2 and F was defined after Eq. (|23|). Similar expression was obtained earlier in Ref. 11211 



for a more involved regime of strong spin-orbit coupling. The last term of (B.9i is purely imaginary and, at low 
frequencies (adiabatic condition \cl>„,\ <s: Oq), corresponds to the Berry phase. Substituting EqB.9 into EqB.8 we 
obtain 



where 



/? 



>S(p = const. — -- — h S^ii 
47* 



/3 ySr (Do 



(B.IO) 



^27 n 12-/ 2 (I>2 + w2| 

m>0 m>0 ^ m/ 

pr cj, 



{6^^.„,6Q>,,-,„ + 6<i>y.„,6^y,-„,) 



m>0 "^ '"ni 

At this point we can perform the Gaussian integration over 50. We obtain 



(B.ll) 



oji + 0)2 



N' 



(l-v7)sinh[^] 



+ (l-v7)2(&2 sinh[^(l-y7)] 



(B.12) 



where N' is a 7-independent and <l>o-independent normalization. 

Thus we see that close enough to Stoner instability, i.e., for 7/7* = (1 - v7) -^ 0, we reproduce Eqs. ( 37 1 and (J38 
Indeed, this is so if we are allowed to replace 



(1 - vJ) 



sinh[^(l-v7)] /3^o 



(B.13) 



for all relevant Oq. This approximation works, at least, in the regime of relatively high temperatures, i.e., when 
V" ' ~ 7 <s r <; 7* . In this case the integral ( 3 8 1 is dominated by Oq ~ ■/* and we obtain jSOq ( 1 - v7) ~ y67 « 1 . 
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Appendix B.3. Restoring the Goldstone mode in Cartesian coordinates 

Having already reproduced the results of the main text we would like to improve our understanding as well as the 
precision of the calculation by analyzing Eq. |B.10 and Eq. B. 1 1 a bit closer. We rewrite Eq. B. 1 1 as 



= - Z 



2J 



6^.„,6<b- 



:„ni^^^-m,z 



m>0 ^ ' 



_ V^ 



,tf) 2 '^0(<^l + O 



[6^-,,mS^x,-m + (y<l)y,„,(5<l)^, _„,) 



n>0 



m>0 » 



(5C)^,,„50,, -m - 50,,„,(5<1).,,_,„) . 



(B.14) 



The second term of the RHS of ( B. 14 1 is problematic since it gives a finite "mass" for the transverse fluctuations. Yet, 
from the spherical symmetry we should expect a Goldstone mode and no "mass". However, the origin of this term is 
quite clear It can be nicely rewritten as 

_ A ^ /.- -- , .^ .^ ^ _ ^ 

47* 



- ^ [6<^x,m6^.x.-m + 6%.,,n5^y,-„^ 



mi:!.) 



47* ° 



(B.15) 



Here AOq is the prol ongati on of the zero mode vector Oq due to transverse fluctuations. This term combines with the 
term -/3<^l/(AJ*) of ( B. 10 1 to give -j63)y/(47*), where O^ = O^ + AO^. Thus we "loose" the Goldstone mode because 
we use Cartesian coordinates and not the spherical ones where Oq would be kept constant. This term could be just 
subtracted thus recovering the Goldstone physics and reducing Oq back to (Do- We distinguish two regimes: 

(A) u„ » Oo 



At high frequencies the action ( B. 14 1 reduces (with no subtractions) to 



(B.16) 



In this limit the integration over the Gaussian fluctuations produces a fac tor (J//3)^ per Matsubara frequency, 
which compensates the factor (J3/J)^ per Matsubara frequency in Eq. (B.5 1. Thus it is not necessary to subtract 
anything at high frequencies. 
(B) w„, « Oo 

In this limit, after the subtraction of the second term of the RHS of ( B. 14 1 and retaining only transverse fluctu- 
ations we obtain 



'(50j_ - 



- Z 






^-f 2 

in>0 



m>0 



2O0 



^l 



In the time representation Eq |B.17| reads 



p 



dr 






40, 



/ 



dr 



(m,f + ((?(Dy)^ 



(B.17) 



(B.18) 



The first term in Eq B.18 is equal to the Berry phase term (23 1. The second term is of the second order in 
derivatives and can be rewritten as 

- [dT(it)\ 

t'O J 



r 

4<t>o 



(B.19) 
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where ft is the unit vector in the direction of <E). Thus we reproduce the second order term (27 1. Since the 
obtained action is local in time and is proportional to the time derivatives, the fact that we started from an 
expansion in small fluctuations around a given Oq is no longer important. Indeed, we can always choose the 
direction of <Po close to the "current" 0(t). 

Appendix B.4. Calculation with a proper counter term 

As we have seen we have to subtract the "mass" term only for the adiabatic frequencies. Here we propose a 
particular form of the counter term which satisfies this condition. We rewrite again Eq. B. 1 1 as iS^o 
where 



Qren , Qcoimfer 



Qren 



z 

m>0 



6<i>,,^„6^-,n,z 



£ 

27 

2/0)2+0.2 



((50.i,,„5(5,,-,„ + 50,,„,(J<1>,,_,„) 



y^ 



L0„ 



m>0 



and 



^ 27* 



27* 0)2 + ^2 V > I 



P 



H 



(B.20) 



(B.21) 



The splitting is chosen so that S'i'^'""''' represents the "mass" term at adiabatic frequencies but vanishes at high fre- 



quencies. We subtract the counter term >S^^""'°'^. Then we obtain 



g)"/,..*,, 



"■n 



N 2 



.0.2 - 



m>0 



(^t 



n 

m>{) 



col + (Jy^o)^ 



N 



,sinh[^] 



2 



7vsinh[^] 



(B.22) 



sinh[^] 
Close to Stoner instabilty, i.e., for 7v — > 1, we can approximate the second multiplier (in rectangular brackets) by 



1 and we reproduce the result (38 1 and (39i from there. In this approach no limitation on temperature from below 
appears. 



Appendix C. QD in presence of an external magnetic field B 

Here we provide an alternative derivation of the susceptibility by considering the isolated QD in a weak magnetic 
field. Unlike in the previous calculation, here the normalization factors of the path integral are not important. 
In the presence of a constant magnetic field B the Hamiltonian of the QD can be written as 



H = ^ 60.4 ^fl„^ - 7S2 - g ■ S , 



(C.l) 



Following the same procedure as described in Appendix B.l we derive the effective action for this situation. The 
constant magnetic field adds to the zero frequency component of the exchange field <Po and, eff'ectively, the new field 
A = (<S>(, + B) plays the role of Oq in the calculation. We obtain 



»So,B - 



m 



2^ trln^[(G^„„ j5„,,„, - 5<D(„,_„,, ' SJ - ^ - 2^ ^ 



m#0 



(C.2) 



where in G^\^ = /e„ - (€„ - p) - AS , is the single particle Green's function of the electrons in presence of B and 
A - JOq + b2 + 2O0B cos 6. Here is the angle between B and <5o. 
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We again expand S(j,^b up to the second order in fluctuations and obtain 

>So,B ~ const. - — + — 2j T7 Zj 2 i^^ ' ^^^Ai^^b ■ S)J . (C.3) 



Foflowing the same steps as described in Appendix B.2 and Appendix B.3 we obtain the partition function 



Z^n(j] ft/fJoOo r<i6lsin6iexp 



The susceptibility is now given by 



47 



exp 



/3vA- 



sinh[f] 



(C.4) 



X = - 



T d^Z\ 



(C.5) 



Performing the appropriate saddle point approximation (up to the quadratic fluctuations around the saddle point) we 



obtain the susceptibility given in Eq. 39 
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